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Abstract 



Q-i, The viability of the lightest neutralino as a dark matter candidate in the 

Next-to-Minimal Super symmetric Standard Model is analysed. We carry out a 
thorough analysis of the parameter space, taking into account accelerator con- 
straints as well as bounds on low-energy observables, such as the muon anomalous 
magnetic moment and rare K and B meson decays. The neutralino relic den- 
sity is also evaluated and consistency with present bounds imposed. Finally, the 
neutralino direct detection cross section is calculated in the allowed regions of 
the parameter space and compared to the sensitivities of present and projected 
dark matter experiments. Regions of the parameter space are found where ex- 
perimental constraints are fulfilled, the lightest neutralino has the correct relic 
abundance and its detection cross section is within the reach of dark matter 
detectors. This is possible in the presence of very light singlet-like Higgses and 
when the neutralino is either light enough so that some annihilation channels are 
kinematically forbidden, or has a large singlino component. 



PACS: 12.60. Jv, 95.35. +d 



1 Introduction 



One of the most interesting candidates for dark matter, within the class of Weakly 
Interactive Massive Particles (WIMPs), is the lightest neutralino (x?), which arises 
in the context of R-parity conserving supersymmetric (SUSY) models [1]. Although 
neutralino dark matter has been extensively examined in the framework of the Minimal 
Supersymmetric Standard Model (MSSM) [2], this is not the case of the Next-to- 
Minimal Supersymmetric Standard Model (NMSSM) [3]. 

The NMSSM is an extension of the MSSM by a singlet superfield S, which provides 
an elegant solution to the so-called /i problem of the MSSM [4]. In this case, an 
effective \x term is dynamically generated when the scalar component of S acquires a 
vacuum expectation value (VEV) of order of the SUSY breaking scale. In addition, and 
when compared to the MSSM, the NMSSM renders less severe the "little fine tuning 
problem" of the MSSM [5], or equivalently, the non-observation of a neutral CP-even 
Higgs boson at LEP II. The presence of additional fields, namely an extra CP-even 
and CP-odd neutral Higgs bosons, as well as a fifth neutralino, leads to a richer and 
more complex phenomenology [6-10]. In particular, a very light neutralino may be 
present [7]. The upper bound on the mass of the lightest Higgs state is larger than in 
the MSSM [8]. Moreover, a very light Higgs boson is not experimentally excluded. All 
these properties also contribute to the emergence of dark matter scenarios that can be 
very different from those encountered in the MSSM, both regarding the relic density 
and the prospects for direct detection. 

The NMSSM must comply with a large number of experimental constraints. In 
addition to LEP II and Tevatron limits on the spectrum, one should also take into 
account SUSY contributions to low-energy observables. The most stringent bounds 
arise from kaon and B decays, as well as from the muon anomalous magnetic moment, 
susy_ N a t U rally, in order to be a viable dark matter candidate, the NMSSM lightest 
neutralino must also satisfy the present astrophysical bounds on the relic abundance. 
It is important to notice that, in addition to the channels already present in the MSSM, 
the NMSSM offers new kinematically viable annihilation channels and potential reso- 
nances [11-13], due to the presence of light scalar and pseudoscalar Higgses. 

The direct detection of neutralino dark matter has been discussed in the framework 
of the NMSSM in Refs. [13-16]. In particular, and as pointed out in [16], the exchange 
of very light Higgses can lead to large direct detection cross sections, within the reach 
of the present generation of dark matter detectors. There are currently a large number 
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of facilities devoted to the direct detection of WIMPs through their elastic scattering 
with nuclei. Among them, it is important to refer to the DAMA collaboration [17, 18] 
(which has already reported data favouring the existence of a WIMP signal 1 ), CDMS 
Soudan [20], EDELWEISS [21] and ZEPLIN I [22], as well as the upcoming detectors 
working with 1 tonne of Ge/Xe [23]. 

The purpose of the present work is to extend the analysis of Ref. [16]. In addition to 
upgrading the experimental constraints already present in the previous work, we now 
include the computation of the SUSY contribution to the muon anomalous magnetic 
moment, and bounds from K- and 5-meson decays. Moreover, we incorporate the 
evaluation of the relic abundance of the lightest neutralino in the NMSSM. 

Our work is organised as follows. We briefly address the most relevant aspects of 
the NMSSM in Section 2, paying especial attention to the minimisation of the Higgs 
potential and to the new features of the Higgs and neutralino sectors. We also comment 
on the various constraints on the low-energy observables that will be included in the 
analysis. In Section 3 we introduce the NMSSM lightest neutralino as a dark matter 
candidate and discuss the potential implications regarding its detection cross section 
and relic abundance. Our results are presented in Section 4, where we first analyse the 
effect of the experimental constraints on low-energy observables, a^ USY and b — > 57, 
on the parameter space (Sec. 4.1). We then include the constraint on the neutralino 
relic density (Sec. 4.2). Finally, taking all these bounds into account, we evaluate 
the neutralino detection cross section on the allowed regions of the parameter space 
(Sec. 4.3), comparing the results with the sensitivity of dark matter detectors. We 
summarise our conclusions in Section 5. 

2 The low-energy NMSSM 

In this section we briefly review some of the most relevant aspects of the NMSSM, and 
summarise the main constraints considered when analysing the parameter space. 

1 Notice that the DAMA result has not been confirmed by the other experiments. For attempts to 
show that DAMA and these experiments might not be in conflict, see Ref. [19]. 
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2.1 Overview of the model 

With the addition of a gauge singlet superfield, the MSSM superpotential is modified 
as follows: 

W = (Y u H' 2 Q l u + Y d H{ Q j d + Y e H{ V e) - e tj \ S H{H 3 2 + \kS\ (2. 1) 

where Hf = (if?, iff), Hi = (H+,H$), i,j are SU(2) indices, and e 12 = 1. In 
this model, the usual MSSM bilinear \i term is absent from the superpotential, and 
only dimensionless trilinear couplings are present in W. However, when the scalar 
component of S acquires a VEV, an effective interaction jiHiH 2 is generated, with 
/j, = X(S). Likewise, the soft SUSY breaking terms are accordingly modified, so to 
include new soft-breaking masses for the singlet (m|), and additional trilinear couplings 
(Ax and A K ). After electroweak (EW) symmetry is spontaneously broken, the neutral 
Higgs scalars develop the following VEVs: 

(H° 1 )=v 1 , (H° 2 )=v 2 , (S) = s. (2.2) 

After computing the tree-level scalar potential 2 , V^Jjj^, one must ensure the presence 
of a true minimum. The minimisation can be conveniently separated into two steps. 
First, one imposes a minimum of V^^rai with respect to the phases (signs) of the 
VEVs, while in a later stage one derives the conditions regarding the modulus of the 
VEVs. The first step immediately allows to directly exclude combinations of signs 
for the parameters. In fact, and working under the convention where A and tan/3 are 
positive, and k, /i(= As), A\ and A K can have both signs, one can analytically show 
that there are only four distinct combinations of signs for k, A\, A k , and ji that ensure 
the presence of a minimum of V^f^i [16]: 

(i) k > 0, sign(s) = sign(A A ) = -sign(A K ), 

which always leads to a minimum with respect to the phases of the VEVs. 

(ii) k > 0, sign(s) = -sign(A A ) = -siga(A K ), 

with \A K \ > 3\v iV2\Ax\/(— \sA\\ + k|s 2 |), where the denominator has to be posi- 
tive. 

(iii) k > 0, sign(s) = sign(A A ) = sign(A K ), 

with \A K \ < 3AwiV 2 |-4a|/(|sA a | + k\s 2 \). 

2 For details about the Lagrangian and the neutral Higgs potential, V^Jj^j, of the model see 
Ref. [16]. 
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(iv) k < 0, sign(s) = sign(A A ) = sign(A K ), 

with \A K \ > 3Xv iv 2 \A\\/(\sA\\ — k|s 2 |), where the denominator has to be positive. 

The additional conditions regarding the minimisation of the potential with respect to 
the Higgs VEVs can be derived, and these allow to re-express the soft breaking Higgs 
masses in terms of A, k, A\, A K , v±, v 2 and s: 

m 2 Hi = - A 2 (s 2 + v 2 sin 2 (3) - \m\ cos 2(3 + As tan (3 (ks + A x ) , 
m 2 H2 = - A 2 (s 2 + v 2 cos 2 (3) + -Aff cos 2(3 + As cot (3 (ks + A x ) , 

m 2 s = - A V - 2k 2 s 2 + Xkv 2 sin 2(3 + XAxV sin 2(3 - kA k s , (2.3) 

where v 2 = v \ + v\ = 2My V /g'i and tan (3 = v 2 /vi. 

The neutral Higgs spectrum in the NMSSM includes three scalars, and two pseu- 
doscalar states, whose mass matrices have been discussed in [16]. For our present 
analysis, let us just recall that in either case interaction and physical eigenstates can 
be related as 

h° a = S ab H° b (a,6=1...3), a° = P^P° = 1,2), (2.4) 

where S (P) is the unitary matrix that diagonalises the 3x3 (2x2) scalar (pseudoscalar) 
Higgs mass matrix. In both sectors, we order the physical eigenstates as m h o < m h o < 
m h o and m a o < m a o . The singlet component of the lightest CP-even Higgs is therefore 
given by S 13 . 

A final comment concerns the neutralino sector, now extended to include the 
singlino component of the new chiral superfield, S. In the weak interaction basis defined 
by ^ oT = (j3° = W® = —iX 3 , H®, H®, S*j the neutralino mass matrix reads 

/ Mi -M z sin W cos (3 M z sin W sin (3 \ 

M 2 M z cos 9 W cos [3 -M z cos6 w sin/3 

— M z sin 6*^ cos /3 M z cos 9w cos (3 —As — \v 2 

Mz sin 0\v sin (3 —Mz cos 9w sin (3 —As — Xv\ 

y -Xv 2 -Xvi 2ks J 

(2.5) 

The above matrix can be diagonalised by means of a unitary matrix N, 

^Aiy-oN^ 1 = diag(m^o, m^o, m^o, m^o, m^o) , (2.6) 

where m^o is the lightest neutralino mass. Under the above assumptions, the lightest 
neutralino can be expressed as the combination 

^> = Nn B° + N 12 W% + N lz Hl + N 14 H° + N 15 S . (2.7) 
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In the following, neutralinos with Nf 3 + Nf A > 0.9, or Nf 5 > 0.9, will be referred to as 
Higgsino- or singlino-like, respectively. 

2.2 Constraints on the parameter space 

In addition to ensuring the presence of a minimum of the potential, other constraints, 
both theoretical and experimental, must be imposed on the parameter space generated 
by the low-energy NMSSM degrees of freedom 

A, k, tan/3, fi, A x , A K . (2.8) 

The soft supersymmetry-breaking terms, namely gaugino masses, Mi.^,3, scalar masses, 
m Q,L,u,D,E, and trilinear parameters, Aq tLtUjDjE , are also taken as free parameters and 
specified at low scale. 

A comprehensive analysis of the low-energy NMSSM phenomenology can be ob- 
tained using the nmhdecay 2.0 code [24]. After minimising the scalar potential, thus 
dismissing the presence of tachyons and/or false minima, the Higgs boson masses are 
computed, including 1- and 2-loop radiative corrections. Squark and slepton masses 
are also calculated, as well as the corresponding mixing angles for the third genera- 
tion. Chargino and neutralino masses and mixings are evaluated and all the relevant 
couplings are derived. 

Even though the general analysis is performed at low-energy, a further theoretical 
constraint can be derived, namely the absence of a Landau pole for A, k, Y t and Y& 
below the GUT scale. Including logarithmic one-loop corrections to A and k, the latter 
constraint translates into A < 0.75, \k\ < 0.65, with 1.7 < tan/3 < 54. 

On the experimental side, nmhdecay 2.0 includes accelerator (LEP and Teva- 
tron) constraints, 5-meson decays, and dark matter relic density through a link to 
micrOMEGAS [25]. In particular, direct bounds on the masses of the charged parti- 
cles {H ± , x ± , Q, and on the gluino mass are taken into account [26,27]. Excessive 
contributions to the invisible decay width of the Z boson [28,29], as those potentially 
arising from Z — > x?X? an d Z — > h°a°, are also excluded from the parameter space. 
Finally, in the neutral Higgs sector, one checks the constraints on the production rates 
for all the CP-even states h° and CP-odd states a , in all the channels studied at 
LEP [30]: e + e~ — > h°Z, independent of the h° decay mode (IHDM); e + e~ — > h°Z, 
dependent on the h° decay mode (DHDM), with the Higgs decaying via h° — > bb, 
h° — > t + t _ , h° — > 2 jets h° — > 77 and h° — * invisible; associated production modes 
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(APM), e+e~ -> /i°a°, with h°a° -> 46's, /i°a° -> 4r's and /i°a° -> a°a°a° -> 66's. 
In addition to the latter, one also takes into account the possible two body decays of 
all CP-even, CP-odd and charged Higgs bosons into squarks and sleptons, as well as 
radiatively induced decays of neutral Higgs bosons into two photons and two gluons. 

Regarding I?-meson decays, and although the nmhdecay 2.0 code already contains 
a rough estimate of the b — > 57 decay branching ratio (evaluated at the leading order 
in QCD), we include in our code a more precise computation of the b — > 57 decay in 
the NMSSM [31], taking into account next-to-leading order (NLO) contributions [32, 
33], following the results of [34]. However, we only include leading order (LO) SUSY 
contributions to the Wilson coefficients at the Mw scale 3 . The calculation within the 
context of the MSSM at LO and NLO can be found in [35] and [36], respectively. The 
most recent experimental world average for the branching ratio (BR) reported by the 
Heavy Flavour Averaging Group is [28, 37] 

BR exp (6 -> si) = (3.55 ± 0.27) x 10~ 4 . (2.9) 

On the other hand, the current SM prediction for the branching ratio is [38] 

BR SM (6 -> s 7 ) = (3.73 ± 0.30) x 10~ 4 , (2.10) 

where the charm-loop contribution has been included [33]. We have estimated the 
theoretical error that results from varying the scales in the b — > 57 calculation within 
the NMSSM, following the method described in [34]. We add to this the experimental 
error in quadrature. This procedure is performed at every point of the parameter space, 
typically leading to an error of about 10% of the total BR(6 — > S7) value. Consistency 
at 2a with the experimental central value of Eq. (2.9) is then demanded. 

We have also included in our code other constraints coming from the contribution 
of a light pseudoscalar a in NMSSM to the rare B- and i^-meson decays [31]. When 
the pseudoscalar is very light it could be produced in meson decays and significantly 
affect the rates for K — K and B — B mixing and other SM decays. In particular, our 
code takes into account the constraints from the pseudo-scalar indirect contributions 
to K - K and B - B mixing, B -> B -> X s /! + /i~, B~ -> K~vv, B -> K° S X°, 

3 No charm- loop contributions were included in the analysis of Ref. [34] , giving a SM central value 
of BR SM (6 — > sj) = 3.293 x 10 -4 . This result is obtained by extrapolating the value of the branching 
ratio evaluated at S = 0.9 and /if, = me,, where S parameterises the photon energy cut _E 7 > (1 — 5)rrib/2 
and Hf, is the renormalisation scale. The corresponding new physics contribution has been implemented 
in our code by using the parameterisation of [34] evaluated at /x& = m b and S — 0.9. 
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and by the direct production, at large tan/3, via b — > sa°, — > Xa°, and 5 — > 7ra° 
decays. 

Finally, in our analysis we will also include the constraints coming from the SUSY 
contributions to the muon anomalous magnetic moment, = (g^ — 2) [39]. Taking into 
account the most recent theoretical predictions for this quantity within the SM [40-42] 
and the measured experimental value [43] , the observed excess in a^ xp constrains a pos- 
sible supersymmetric contribution to be a^ USY = (27.6 ± 8) x 10~ 10 , where theoretical 
and experimental errors have been combined in quadrature. 

The evaluation of a^ USY in the NMSSM has been included in our analysis, and 
consistency at the 2a level imposed. Thus those regions of the parameter space not 
fulfilling 11.6 x 1(T 10 < a£ USY < 43.6 x KT 10 will be considered disfavoured. 



3 Dark matter in the NMSSM 

The new features of the NMSSM have an impact on the properties of the lightest 
neutralino as a dark matter candidate, affecting both its direct detection and relic 
abundance. 

The computation of the spin-independent part of the neutralino-nucleon cross sec- 
tion was discussed in detail in [16]. It was pointed out there that the existence of a 
fifth neutralino state, together with the presence of new terms in the Higgs-neutralino- 
neutralino interaction (which are proportional to A and k) , trigger new contributions to 
the spin- independent part of the neutralino-nucleon cross section, cr^o_ p . On the one 
hand, although the term associated with the s-channel squark exchange is formally 
identical to the MSSM case, it can be significantly reduced if the lightest neutralino 
has a major singlino composition. On the other hand, and more importantly, the dom- 
inant contribution to cr^o_ p , associated to the exchange of CP-even Higgs bosons on 
the t-channel can be largely enhanced when these are very light. In the NMSSM, the 
lightest CP-even Higgs can escape detection if its singlet composition is large. For 
instance, this makes possible the presence of scenarios with m h o < 70 GeV, thus con- 
siderably enhancing the neutralino-nucleon interaction. Consequently, large detection 
cross sections can be obtained, even within the reach of the present generation of dark 
matter detectors. 

However, in order to be a good dark matter candidate, the lightest NMSSM neu- 
tralino must also comply with the increasingly stringent bounds on its relic density. 
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Astrophysical constraints, stemming from the analysis of galactic rotation curves [44], 
clusters of galaxies and large scale flows [45] , suggest the following range for the WIMP 
relic abundance 

0.1<^ 2 <0.3, (3.11) 

which can be further reduced to 

0.095 < Vlh 2 < 0.112, (3.12) 

taking into account the recent three years data from the WMAP satellite [46]. 

Compared to what occurs in the MSSM, one would expect several alterations re- 
garding the dominant processes. As discussed in [13], and as mentioned above regarding 
the direct detection cross section, the differences can be present at distinct levels. First, 
and given the presence of a fifth neutralino (singlino), the composition of the annihi- 
lating WIMPs can be significantly different from that of the MSSM in wide regions of 
the parameter space. Having the possibility of a singlino-like lightest supersymmet- 
ric particle (LSP), associated with the presence of new couplings in the interaction 
Lagrangian, in turn favours the coupling of the WIMPs to a singlet-like Higgs, whose 
mass can be substantially lighter that in the MSSM, given the more relaxed experimen- 
tal constraints. Regarding the channels through which neutralino annihilation occurs, 
in the NMSSM we have new open channels, essentially due to the existence of light 
Higgs states. In summary, the presence of additional Higgs states (scalar and pseu- 
doscalar) favours annihilation via s-channel resonances. On the other hand, having 
light hi and states that are experimentally viable means that new channels with 
annihilation into Z h\, h\ h®, h\ a\ and a? a\ (either via s-channel Z, h®, a° exchange or 
t-channel neutralino exchange) can provide important contributions the annihilation 
and co- annihilation cross-sections [12]. 

Noticing that important annihilation channels (s-channel) are related to the t- 
channel processes responsible for the most relevant contributions to <J^,o_ p , one should 
expect a strong interplay between a viable relic density, and promising prospects for 
the direct detection of the NMSSM dark matter candidate. In fact, there should be 
regions of the parameter space which provide new and interesting scenarios 4 . 

4 As concluded in [13], it might even be possible to reconcile a very light neutralino with the 
experimental observations from DAMA, CDMS II, and WMAP. 
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4 Results and discussion 



In this section, we study the viability of lightest NMSSM neutralino as a good dark 
matter candidate. Motivated by the results obtained in [16], we focus on regions of 
the low-energy NMSSM parameter space where large direct detection cross sections are 
likely to be obtained. Building upon the previous analysis, we apply the new constraints 
(improved comparison with LEP and Tevatron data), K- and S-meson decays, a^ USY , 
and compatibility with the WIMP relic density. Finally, we discuss the prospects of 
the experimentally viable regions regarding direct detection of dark matter. 

Let us just recall that the free parameters of the model, associated with the Higgs 
and neutralino sectors of the theory, are 5 

A, k, n(=\s), tan/3, A x , A K , M 1 , M 2 , M 3 . (4.13) 

We assume that the gaugino mass parameters mimic, at low-energy, the values of a 
hypothetical GUT unification = M 1 = 

It should be emphasised that in nmhdecay 2.0 some of the input parameters are 
specified at a different scale than in the former version nmhdecay 1.1 [9], which was 
used in the previous analysis [16]. Although the difference between the values of A 
(or k) at the EW and the SUSY scales (~ 1 TeV) is very small, there is a substantial 
change in the value of the trilinear coupling A\. These variations are induced by the 
top trilinear coupling A top , and are approximately given by A^ USY m Af w + Af™. 
Therefore, one needs to take this shift into account when comparing the present results 
with those of [16]. 

Motivated by the results of [16] regarding the prospects for direct detection of dark 
matter, we will be interested in a regime of low tan /3, as well as in values of in the 
range 110 GeV < fi < 200 GeV (the lower limit ensuring that in most cases one can 
safely avoid the LEP bound on the lightest chargino mass). Likewise, the following 
intervals for the trilinear couplings will be taken: —800 GeV < A\ < 800 GeV, 
and —300 GeV < A K < 300 GeV (the optimal ranges will typically correspond to 
\A\\ ~ 400 GeV and \A K \ < 200 GeV, working in a small tan/5 regime). 

Slepton and squark masses, as well as the corresponding trilinear parameters, do 
not significantly affect the neutralino detection properties, other than through the 

5 Although the soft gluino mass, M 3 , is not directly related to the computation of the 
Higgs/neutralino masses and mixings, it plays a relevant role in contributing to the radiative cor- 
rections to the Higgs boson masses. 
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radiative corrections to the Higgs masses. However, low-energy observables are very 
sensitive to their specific values. In the following section we will see, for instance, how 
the experimental constraint on a^ USY favours light sleptons. 

As already done in Ref. [16], we divide the scan of the low-energy NMSSM param- 
eter space following the results of the minimisation with respect to the VEV phases, 
separately discussing each of the cases (i)-(iv) (see Section 2.1). 



4.1 NMSSM parameter space: updated constraints 

We first discuss the new constraints on the parameter space arising from the improved 
analysis on the Higgs sector, the muon anomalous magnetic moment, and K- and 
5-meson decays. 

Among the new features implemented in nmhdecay 2.0, one finds additional ra- 
diative corrections to the Higgs boson masses, including corrections of order g 2 Y 2 b to 
the CP-even Higgs boson mass (induced by stop/sbottom D-term couplings). Regard- 
ing the logarithmic one-loop corrections of the order g 2 , these are now dependent on 
the different masses of squarks/sleptons of distinct generations. Moreover, the correc- 
tions to fourth order in A and k are also taken into account. The computation of the 
sparticle spectrum is also complete in the new version, and all squark and gluino data 
is confronted with the constraints from both Tevatron and LEP. With respect to the 
results obtained in the previous analysis [16], the latter improvements only translate 
into slight changes in the exclusion regions. 

Concerning the evaluation of the supersymmetric contributions to the muon anoma- 
lous magnetic moment, the relevant processes comprise neutralino-sneutrino as well as 
chargino-smuon loops. The only change with respect to the MSSM is due to the 
fifth neutralino state and the corresponding modified neutralino-lepton-slepton cou- 
pling. Since we are interested in cases with very low tan/5, the contributions from 
neutralino and chargino loops are of similar magnitude, and very small. For example, 
with tan (3 = 3 and slepton mass parameters above rriE,L ~ 1 TeV one typically obtains 
C V ~ 10~ n , which is disfavoured. In order to obtain compatibility with the experi- 
mental result an increase in the value of tan f3 is welcome, but this would then lead to 
regions of the parameter space where, from the dark matter point of view, the NMSSM 
resembles the MSSM. The other possibility is decreasing the slepton (and gaugino) 
masses. Furthermore, large values of the slepton trilinear couplings are needed in or- 
der to increase the LR mixing in the smuon mass matrix. The choice [/,Ae < is 
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optimal, since it makes the neutralino contribution positive and large. For example, 
with tan/5 ~ 5, A E = -2500 GeV and m EjL < 200 GeV, one obtains a^ USY > 10~ 9 
for Mi < 215 GeV. The relevance of these changes is illustrated in Fig. 1, where the 
numerical results for a^ USY are plotted versus the bino mass, Mi, for different combi- 
nations of slepton mass and trilinear couplings. For each case we have also varied A\ 
and A K over a wide range and scanned the whole (A, k) plane, which as evidenced in 
the figure has virtually no effect on the resulting a^ USY . We have also included the 
various LEP and Tevatron constraints. For the rest of our analysis, we will assume 
m E ,L = 150 GeV and A E = —2500 GeV. Also, and unless otherwise stated, we will set 
the bino mass to Mi = 160 GeV, which, according to Fig. 1, leads to a sufficiently large 
a susY_ detection properties of the neutralino are in general quite insensitive to 
changes in the slepton sector. Notice however that if one does not wish to impose the 
bound on the muon anomalous magnetic moment, heavy sleptons (equal to squarks) 
can be taken which would not affect the dark matter predictions. 

Regarding the bounds arising from K- and £>-meson physics, the most important 
role is played by the b — > S7 decay, which can in principle exclude important regions 
of the parameter space. Concerning the other K- and £>-meson processes discussed 
in Section 2.2, we have verified that throughout the investigated NMSSM parameter 
space they are always in good agreement with experiment, so that we will make no 
further reference to the latter in the following discussion of the numerical results. 

In the present analysis we will not take into account any source of flavour violation 
other than the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Moreover, we will be 
systematically considering large values for the gluino mass (above 1 TeV). Under the 
latter assumptions, the most important contributions to BR(5 ^57) arise in general 
from charged Higgs and chargino mediated diagrams [35]. 

On the one hand, when the dominant contributions are those stemming from 
charged Higgs exchange, the results for BR(5 ^57) closely follow the behaviour 
of the charged Higgs mass, which in the NMSSM is given by 

2 2 /i 2 K 2 2 2 [L A\ 2 t A -\ A\ 

mH± = ^W)x - v x + Mm + w ' (414) 

From the above, we expect that smaller values of BR(6 — > s 7) should be obtained for 
large m 2 H ±, and therefore when k/X is sizable (for positive values of k) or for small 
k/X (if k < 0). In general, smaller values of the BR(6 — > 57) will be also associated 
to larger values of the product fiA\. Furthermore, the leading term of the Wilson 
coefficient associated to the charged Higgs varies as tan -2 /? [35]. As a consequence, 
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Figure 1: Supersymmetric contribution to the anomalous magnetic moment of the muon 
as a function of the bino mass, M lr for tan/3 = 5, fx = 150 GeV, trilinear couplings in 
the range -800 < A x < 800 GeV, -300 < A K < 300 GeV. From bottom to top, the 
different bands correspond to the following values of the slepton mass and lepton trilinear 
terms, m L:E = 1 TeV with A E = 1 TeV, m L:E = 150 GeV with A E = 1 TeV, and 
m L:E = 150 GeV with A E = —2.5 TeV. A full scan on the (A, k) plane has been performed 
for each case, including LEP and Tevatron experimental constraints. The horizontal solid 
line indicates the lower bound of the allowed 2a interval. 

one expects a decrease of this contribution as tan (3 increases. On the other hand, in a 
regime of \i < M 2 , the lightest chargino is Higgsino-dominated, so that its mass is also 
quite small (ra-± ~ ji). Thus, the chargino contributions (which are opposite in sign to 
those of the charged Higgs) are also expected to play a relevant role, although, in the 
cases analysed in this paper (tan/? < 10), they are not dominant. Gluino contributions 
are also very small, given the little flavour violation in the down squark sector, and the 
sizable values of M 3 . Likewise, neutralino exchange contributions are almost negligible. 
We thus find that, in general, the NMSSM contribution to BR(6 — > s^y) at low tan/3 is 
large and mostly arising from charged Higgs loops. This leads to stringent constraints 
on the parameter space. 

Let us study the effect of the experimental bound on BR(6 — > 57), together with 
the updated accelerator constraints on the NMSSM parameter space. After this first 
survey we will no longer separately address the K- and _B-meson constraints (and 
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a M ) from those arising from LEP/Tevatron data. Henceforth, experimentally allowed 
regions will be those that not only comply with the latter data, but that also exhibit 
BR(6 — > S7) within 2cr from its central experimental value. As mentioned before, 
in order to satisfy the constraint on the muon anomalous magnetic moment, we take 
Mi = 160 GeV in the following subsection, for which, in the case with tan/3 = 5, 
susy w 14 x 1Q -9 ( gee pig. 1). 



4.1.1 fiA K < and /iA x > (re > 0) 

As discussed in [16], this is one of the most interesting areas of the parameter space, 
since although sizable regions are excluded due to the occurrence of tachyons in the 
CP-even Higgs sector (namely for larger values of \A K \), the possibility of having ex- 
perimentally viable light scalar Higgs leads to potentially large values for a^o_ p . 

As an example, we represent on the left-hand side of Fig. 2 the (A, re) parameter 
space for an example with tan (3 = 3, A\ = 200 GeV, A K = —200 GeV and /i = 
130 GeV. The tachyonic region in the CP-even Higgs sector is depicted, as well as the 
region excluded due to the presence of false minima of the potential. An important part 
of the theoretically allowed region is also ruled out due to conflict with LEP and/or 
Tevatron data. This owes to the fact that the doublet component of the lightest scalar 
Higgs is very large and gives rise to excessive Higgs production rates, in particular, 
e + e~ -> h°Z, IDHM and DHDM (h° -> bb and to a lesser extent, h° -> 2 jets). Once 
all these bounds are applied, a small area on the right of that experimentally excluded 
survives, remarkably exhibiting very light Higgses and neutralinos (associated with 
a singlet /singlino component above 90%) and therefore clearly characteristic of the 
NMSSM. We recall that these are the regions where one expects to find large theoretical 
predictions for a^o_ p . 

On the right-hand side of Fig. 2 we superimpose the results for the BR(b — > S7) on 
the (A, re) plane. As discussed in the previous section, the resulting branching ratio is 
typically large, BR(6 ^ 57) > 3.5 x 10~ 4 , and increases to as much as ~ 5 x 10~ 4 in 
regions with small re/A, where the charged Higgs mass is smaller. Notice therefore that 
the regions of the (A, re) plane associated with larger values of BR(6 — > S7) typically 
correspond to those where the largest predictions for a^o_ p are found. In this example, 
only a small triangular region with A < 0.05, for re < 0.7, is within a la deviation from 
the experimental bound of Eq. (2.9) and A < 0.35 is needed in order to be within 2a 
of that result. In the plot we also indicate with dot-dashed lines the different values of 
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Figure 2: Effect of the experimental constraints on the (A, k) plane for an example with 
tan/? = 3, A\ = 200 GeV, A K = -200 GeV and fi = 130 GeV. In both cases, the 
gridded area is excluded due to the appearance of tachyons, while the vertically ruled area 
corresponds to the occurrence of unphysical minima. The oblique ruled area is associated 
with points that do not satisfy the LEP and/or Tevatron constraints or where (at least) 
the LEP bound on direct neutralino production is violated. The region above the thick 
black line is disfavoured due to the occurrence of a Landau pole below the GUT scale. 
On the left plot, from top to bottom, solid lines indicate different values of the lightest 
Higgs mass, m h o = 114, 75, 25 GeV. Dashed lines separate the regions where the lightest 
scalar Higgs has a singlet composition given by S^ 3 = 0.1, 0.9 (from top to bottom). 
Finally in the area below the dotted line, the lightest neutralino has a singlino composition 
greater than Nf 5 = 0.1. On the right, grey areas represent the theoretical predictions for 
BR(b — > 57). From left to right, la (dark), 2a (medium) and excluded (light) regions 
are shown. Dot-dashed lines stand for the different values of the charged Higgs mass, 
m H ± = 1000, 500, 450 GeV (from left to right). 

the charged Higgs mass, thus illustrating the correlation between its decrease and the 
increase in BR(6 — > 57). 

The effect of the various experimental constraints on the parameter space is very 
sensitive to variations in the input parameters. We will now investigate how changes in 
tan/3 and A\ affect the resulting BR(6 — > S7). On the one hand, as already mentioned, 
increasing the value of tan/3 leads to a reduction of the charged Higgs contribution. 
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Figure 3: (A, k) parameter space for tan/? = 5, A K — —200 GeV, and fi = 130 GeV. 
On the left we take A x = 200 GeV, while on the right we consider A x = 400 GeV. Line 
and colour code follow the conventions of Fig. 2. In this case, the single dot-dashed line 
corresponds to m H ± = 1000 GeV. 

Since in our case this is the leading contribution to BR(6 — > S7), an enhancement 
in tan (3 enlarges the regions of the parameter space which are consistent with the 
experimental constraint. This is illustrated on the left-hand side of Fig. 3 with the 
same example of Fig. 2, but now taking tan (3 = 5. The resulting charged Higgses 
are heavier (m H ± > 500 GeV) and as a consequence the entire (A, k) plane fufils the 
experimental constraint on BR(6 — > 57). Notice that LEP and Tevatron constraints 
are also modified. On the other hand, an increase in the trilinear term A\ also leads 
to heavier charged Higgses as seen in Eq. (4.14). Therefore, this can induce a further 
decrease in BR(6 — > 57). An example of this is shown on the right-hand side of Fig. 3, 
where in addition to tan/3 = 5, A\ = 400 GeV has been used. Again, the whole (A, k) 
plane is allowed due to the increase in m H ±. 

Let us finally comment on the possibility of changing the signs of /i, A\, and A K , 
while keeping /j,A\ > and fiA K < 0. Although the Higgs potential is invariant 
under this change, the same does not occur for the chargino and neutralino sectors, 
so that both these spectra, as well as the experimental constraints are likely to be 
modified. As an illustrative example, we present in Fig. 4 the same case as in the 
left-hand side of Fig. 3 but with the opposite signs for /i, A\, and A K . There are some 
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Figure 4: (A, k) parameter space for tan/3 = 5, A x = —200 GeV, A K = 200 GeV and 
fj, = —130 GeV. Line and colour code following the conventions of Fig. 2. 

important alterations to the areas excluded by unphysical minima and experimental 
constraints, both of which are now more extensive. Finally, notice that the BR(6 — > s 7) 
now excludes a larger area of the (A, k) plane, thereby disfavouring those areas which 
potentially lead to larger neutralino detection cross sections. 

In the light of this analysis, the optimal areas of the parameter space correspond 
to those with /i, A\ > 0, and A K < 0, and where tan/3 and A\ are relatively large. 
In order to keep within the context where NMSSM-like dark matter scenarios can be 
obtained, we will use tan/3 < 5. 

4.1.2 fj,A K < and fj,A x < (k > 0) 

Compared to the previous case, the presence of tachyons gives rise to far stronger con- 
straints. Occurring now in both CP-even and CP-odd Higgs sectors, the non-physical 
(tachyonic) solutions exclude very large areas of the parameter space. Regarding the 
LEP experimental exclusions, these arise from excessive contributions to h° — > bb and 
h° — > 2 jets, and cover an area wider than what had been previously identified in [16] 
(a consequence of the improved computation of the Higgs spectrum). In addition, due 
to the lightness of the charged Higgs bosons, an important region is also excluded due 
to very large BR(6 — > 57), so that the only surviving regions are those associated with 
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Figure 5: (A, k) parameter space for tan/3 = 5, A x = —200 GeV, A K = —200 GeV and 
\i = 130 GeV. Line and colour code following the conventions of Fig. 2. 

A < 0.2. As an example, Fig. 5 displays a case with tan/5 = 5, A\ = —200 GeV, 
A K = —200 GeV and fi = 130 GeV. Let us remark that since the lower-right cor- 
ner of the (A, k) plane is not accessible, one cannot find light neutral Higgs states, so 
that interesting prospects regarding the direct detection of dark matter should not be 
expected. 

Although varying the several parameters results in modifications of the excluded 
areas (LEP/Tevatron, b — > s^y and unphysical minima), in all cases these are sizable. 
Only very reduced regions, corresponding to small values of A survive all the constraints. 
In these areas the singlet component of the lightest Higgs is negligible and the lightest 
neutralino is Higgsino-like, therefore resembling MSSM scenarios. The complementary 
region, with A\, A K > and (j, < 0, leads to even more extensive tachyonic regions, 
and we will not further discuss it. 

4.1.3 fiA K > and jiA\ > (k > 0) 

This combination of signs leads to a parameter space which is plagued with tachyons 
[16], arising from both Higgs sectors. Contrary to what was noticed for the previous 
cases, here the unphysical minima occur for small values of A. The remaining areas in 
the (A, k) plane are also very affected by experimental constraints. As a consequence, 
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Figure 6: On the left (A, k) parameter space for tan (3 = 3, A x = 200 GeV, A K = 50 GeV 
and fj, = 160 GeV. On the right tan/3 = 5, A\ = 400 GeV, A K = 200 GeV, and fj, = 
130 GeV for k < 0. Line and colour code following the conventions of Fig. 2. 

only very reduced areas of the parameter space survive. 

For example, considering the choice A x = 200 GeV, A K = 50 GeV and // = 160 GeV, 
with tan (3 = 3, we observe that once the areas corresponding to the occurrence of 
tachyons are excluded, the small surviving region is still plagued by false minima as 
well as by the violation of several experimental constraints. In particular, lower values 
of k are ruled out due to conflict with LEP (h° — > bb and h° — > 2 jets) and excessive 
contributions to the BR(6 — > 57). This is illustrated on the left-hand side of Fig. 6. 

Although reducing the values of A K and tan (3 enlarges the areas where physical 
minima can be found, the addition of the experimental constraint on BR(6 — > sj) to 
the LEP and Tevatron bounds typically rules out the whole (A, k) plane. Increasing 
tan/5 in order to reduce the contribution to BR(6 — > 57) is also highly disfavoured 
since the tachyonic region becomes more important. As a consequence, no interesting 
implications for neutralino dark matter detection are expected in this case. 
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4.1.4 \iA K > and fiA x > (k < 0) 



The only viable combination of signs associated with negative values of k is also plagued 
by the appearance of tachyons, both in the CP-even and CP-odd Higgs sectors, towards 
the regions with small A [16]. In addition, experimental constraints also exclude large 
portions of the parameter space in the vicinity of the tachyonic regions. As an example, 
let us mention that for the case A\ = 100 GeV, A K = 50 GeV and \x = 130 GeV, with 
tan/5 = 3, all the parameter space associated with physical minima is ruled out, since 
either DHDM constraints (h° — > bb and h° — > 2 jets) are violated or consistency 
with the BR(6 — > s^y) bound is not achieved. Increasing A x , A K , and tan/5 leads 
to a significant improvement. For example, with A\ = 400 GeV, A K = 200 GeV, 
ji = 130 GeV with tan/5 = 5 some allowed areas are found, as shown on the right-hand 
side of Fig. 6. Nevertheless, LEP/Tevatron experimental constraints together with the 
bound on BR(b s 7) rule out those parts of the parameter space where the Higgs is 
light and more singlet-like, and in which 0^0 _ p can be sizable. The remaining allowed 
regions correspond to a rather small area, in which the lightest Higgs is essentially 
doublet-like, while the lightest neutralino exhibits a strong Higgsino dominance. 



4.2 Neutralino relic density 

The next step in our analysis is to take into account the available experimental data 
on the WIMP relic density. In order to be a viable dark matter candidate, the 
lightest NMSSM neutralino must have an abundance within the ranges presented in 
Eqs. (3.11,3.12). Similar to what occurs in the MSSM, this additional constraint fur- 
ther reduces the regions of the low-energy parameter space. Moreover, and as hinted 
before, one expects that fl^o h 2 will in general lie below the experimental ranges. 

A thorough analysis of the relic density of dark matter in the NMSSM has been 
carried out in [12]. It was found that compatibility with the WMAP constraint is 
possible in the regions where the lightest Higgs is dominated by the doublet components 
and the lightest neutralino is a bino-Higgsino mixture. Apart from possible Higgs 
resonances, compatible values of Q^o h 2 are found for two distinct regions: fi ^> M 2 
and /j, > Mi (similar to the MSSM for small to intermediate values of tan/3). Also 
note that, for regions with /i < Mi, Higgsino-singlino neutralinos with masses below 
My/ can give a relic density within the WMAP range (or larger) essentially because the 
annihilation into Z and W gauge bosons is kinematically forbidden. It is also worth 
noticing that a pure bino LSP also offers interesting scenarios, with a remarkable 
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role being played by s-channel Higgs resonances (else f2^o h 2 tends to be above the 
experimental bound). Additional LSP annihilation via scalar or pseudoscalar Higgses 
can also play a relevant role. 

In order to understand the results for f2^o h 2 one needs to take into account the 
variations in the mass and composition of the lightest neutralino in the (A, k) plane, as 
well as in the Higgs sector. In general, the neutralino relic density will be too small in 
those regions of the parameter space where it is Higgsino-like and increases when the 
neutralino becomes more singlino-like. In addition, one should consider the possible 
existence of resonant annihilation (when twice the neutralino mass equals the mass of 
one of the mediating particles in an s-channel) and the kinematic thresholds for the 
various channels (e.g., annihilations into ZZ, WW, Zh®, h®h®, a°a°, and a^h®). 

Since the goal of our present study is to discuss the potential of NMSSM-like sce- 
narios regarding the theoretical predictions for cr^o_ p , in this subsection we focus on 
those examples of the parameter space having large neutralino detection cross section. 
We investigate to which extent the inclusion of the bound of the relic density further 
constrains the parameter space. As pointed out in [16], these scenarios typically occur 
in association to singlet-like h®, with singlino-Higgsino neutralinos. Let us study one 
example in detail. 

We begin by taking M x = 160 GeV, A x = 400 GeV, A K = -200 GeV, and /i = 
130 GeV, with tan /3 — 5, which according to Figs. 1 and 3 is consistent with the bounds 
on a^ USY and BR(6 — > sj), respectively. The results for the neutralino relic density 
are depicted in the (A, k) plane on the left-hand side of Fig. 7. On the experimentally 
allowed area, grey dots stand for points which, in addition to experimental constraints, 
fulfil 0.1 < Q^o h 2 < 0.3, whereas black dots represent points in agreement with the 
WMAP constraint. Notice that in this case, no points are excluded by LEP or Tevatron 
bounds. 

For large values of k and small A (i.e., on the upper left corner of the plots), the 
lightest neutralino is relatively heavy and has a mixed bino-Higgsino composition, since 
we have chosen // ~ Mi. Due to the large Higgsino component, the neutralino relic 
density is very small, and cannot account for the observed amount of dark matter. 
As we move in the (A, k) plane towards smaller values of k and larger values of A, 
the neutralino becomes lighter and has a larger singlino component (in this example 
Nf 5 < 0.35), and as a consequence, fl^o h 2 increases. As the neutralino mass decreases, 
some annihilation channels become kinematically forbidden, such as annihilation into a 
pair of Z or W bosons when m^o < Mz or m^o < Mw, respectively. We have indicated 
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Figure 7: (A, k) parameter space with information about the neutralino relic density. On 
the left, an example with M l = 160 GeV, tan/3 = 5, A\ = 400 GeV, A K = -200 GeV, and 
fj, = 130 GeV. The gridded area is excluded due to the appearance of tachyons, while the 
vertically ruled area corresponds to the occurrence of unphysical minima. The region above 
the thick black line is disfavoured because of the occurrence of a Landau pole below the GUT 
scale. The oblique ruled area is associated to those points that do not satisfy LEP/Tevatron 
and/or BR(6 ^57) constraints, whereas the bound on a^ USY is fulfilled in the whole plane. 
The dark shaded (cyan) area corresponds to points which are experimentally viable, and 
whose relic density complies with the astrophysical bound of Eq. (3.11). Points in black 
are those in agreement with experimental constraints and WMAP bounds (c.f. Eq. (3.12)). 
The dashed red lines indicate the resonances of the lightest neutralino annihilation channels 
through the second lightest CP-even Higgs, 2m^ = m h o. In the region below the red 
dotted line the lightest neutralino mass is larger than the mass of the lightest Higgs. Along 
the red solid lines the neutralino mass is equal to the Z and W mass (from left to right, 
respectively). On the right, the same example is shown, but with A K = and ji = 150 GeV. 
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these two thresholds in the figure with red solid lines. Below these the resulting relic 
density can be large enough to fulfil the WMAP constraint. Variations in the Higgs 
sector also affect the calculation of the neutralino abundance. On the one hand, the 
mass and composition of the lightest Higgs also vary throughout the (A, k) plane. 
Lighter Higgses with a larger singlet composition are obtained for small values of k. In 
our case, m^o < m h o for large k and small A, but eventually, the Higgs becomes lighter 
and new annihilation channels (the most important being X1X1 ~^ h\h\ and Zh\) 
are available for the neutralino, thus decreasing its relic density. The points where 
m^o = m h o are indicated with a dotted red line in the plot. On the other hand, one 
also needs to take into account the existence of rapid neutralino annihilation with the 
second-lightest CP-even Higgs, when 2m^o = m h o, which is responsible for a further 
decrease in fl^o h 2 . This is indicated in the plot with a red dashed line. 

As we can see, in the present example the correct relic density is only obtained 
when either the singlino composition of the neutralino is large enough or when the 
annihilation channels into Z, W, or h® are kinematically forbidden. Interestingly, 
some allowed areas are very close to the tachyonic border. The neutralino-nucleon 
cross section can be very large in these regions, due to the presence of very light 
singlet-like Higgses (in this example Sf 3 ps 0.99). 

The same example, but now with A K = and /i = 150 GeV is shown on the right- 
hand side of Fig. 7. Once more, in order to reproduce the correct h 2 the neutralino 
has to be either sufficiently light so that some annihilations channels are closed or have 
a large singlino component. In this particular case the singlino component of Xi can be 
even larger, with N 2 5 ~ 0.9 in the allowed area with very low k. Notice, however, that 
the region in the vicinity of the tachyonic area is excluded by experimental bounds. 

In order to study the importance of the neutralino composition, we will now consider 
variations in the gaugino masses. To begin with, we increase the bino mass and take 
Mi = 330 GeV, thereby decreasing the bino component of the lightest neutralino. 
Such an increase of the gaugino masses implies a reduced contribution to the muon 
anomalous magnetic moment. We obtain af SY « 7.2 x 10 10 (see Fig. 1), more than 
2a away from the central value and therefore disfavoured. The resulting (A, k) plane is 
represented on the left-hand side of Fig. 8. Since the Higgsino component has increased 
with respect to the previous examples, the resulting relic density for the neutralino in 
the region with large k and small A is even smaller. Once more, in order to have the 
correct Q^o h 2 we need to go to regions of the parameter space where some annihilation 
channels are not kinematically allowed and/or the neutralino is singlino-like. Notice 
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Figure 8: On the left, the same as in Fig. 7 but with Mi = 330 GeV, tan/3 = 5, A x = 
570 GeV, A K = —60 GeV, with /x = 160 GeV. The area to the right of the red dotted line 
has m^o < m h o. Along the red dot-dashed line, 2m^o = Mz, leading to rapid neutralino 
annihilation. On the right, the same with M 1 = 500 GeV, tan/3 = 5, A x = 400 GeV, A K = 
— 150 GeV, with /x = 130 GeV. In this example the red dot-dashed line indicates resonances 
on the neutralino annihilation mediated by the lightest pseudoscalar when 2m^o = m a o. In 
both examples the resulting a^ USY is outside the experimental 2a region. 



also that the neutralino is in general heavier in this example and therefore the lines 
with m^o = Mz and m^o = My/ are shifted to lower values of k. Also, the region with 
m^o < m h Q is modified and now corresponds to the area on the right of the dotted red 
line. Finally, we must take into account the possible resonances along which fl^o h 2 
decreases. In this example, rapid annihilation of neutralinos occurs via CP-even Higgs 
exchange when 2m^o = m h o, which takes place along the two upper red dashed lines. 
There is also a resonance with the Z boson when 2m^o = Mz which occurs along 
the lower red dashed line. It is worth noticing that, once more, a part of the region 
allowed by experimental and astrophysical constraints lies close to the tachyonic area, 
and could have a large cr^o_ p . 

Finally on the right-hand side of Fig. 8 we show another example where the bino 
mass has been further increased to M\ = 500 GeV. The contribution to the muon 
anomalous magnetic moment is also too low. As shown in Fig. 1, a^ USY »4x 10~ 10 , 
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more than 2a away from the experimental value. Due to the further increase of the 
Higgsino component and mass of the lightest neutralino, its relic density is even smaller 
and compatibility with WMAP is only obtained when the neutralino is lighter than 
the lightest Higgs and, at least, the W boson. Notice that in this example there is also 
a resonant annihilation through the lightest CP-odd Higgs when 2 m^o ps m a o , which 
further decreases f2^o h 2 . This constrains the allowed region to small values of A, in 
which the neutralino is mostly singlino, N 2 5 < 0.8. 

To summarise, in these scenarios, neutralinos typically have a very small relic den- 
sity (insufficient to account for the dark matter in the Universe) as a consequence of 
their large Higgsino composition. Only when one moves towards regions of the pa- 
rameter space where the singlino composition is enhanced or the neutralino mass is 
decreased (such that some annihilation channels become forbidden) can the WMAP 
result be reproduced. 

We have not yet addressed the other areas of the parameter space (ii)-(iv). These 
other choices of signs for the different parameters are associated to less favourable sce- 
narios. First, in cases (ii) and (iv) the experimental constraints rule out the regions 
where Xi has a large singlino component. Therefore, the lightest neutralino is in gen- 
eral Higgsino-like throughout all the allowed (A, k) plane and, consequently, its relic 
abundance is much below the favoured values, i.e., Q^o h 2 <C 0.1. In case (iii), where 
large areas are excluded because of the occurrence of tachyons, it is extremely com- 
plicated to find regions that simultaneously fulfil the experimental and astrophysical 
constraints. For this reason, in the following section we will limit our analysis to case 

As a next step of our analysis, we will bring together all the constraints so far 
explored, and after having ensured that we are indeed in the presence of a viable 
NMSSM scenario (namely with the correct relic density), we will investigate to which 
extent the lightest neutralino can be detectable in dark matter experiments. 



4.3 Neutralino direct detection prospects 

After having discussed the new and the improved constraints on the low-energy pa- 
rameter space, we will now address whether or not NMSSM neutralinos with a relic 
density in agreement with current limits are likely to be detected by the present or the 
next generation of dark matter detectors. 
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Figure 9: Scatter plot of the scalar neutralino-nucleon cross section as a function of the 
lightest neutralino mass. On the left, an example with M x = 160 GeV, tan/? = 5, A x = 
400 GeV, A K = —200 GeV, and /i = 130 GeV. All the points represented are in agreement 
with LEP/Tevatron, a^ USY , and BR(6 — > s^y) bounds. Dark gray dots represent points 
which, in addition, fulfil 0.1 < Q^o h 2 < 0.3, whereas black dots are those in agreement 
with the WMAP constraint. The sensitivities of present and projected experiments are also 
depicted, with solid and dashed lines, respectively. On the right we show the same example 
but with fj, = 150 GeV and A K = 0. 

Although in our survey of the low-energy NMSSM parameter space we have scanned 
over all combinations of signs (i)-(iv), as we already mentioned, cases (ii)-(iv) present 
far less interesting situations regarding the neutralino relic density. Even though one 
can find in the latter three cases some challenging situations regarding direct detection 
prospects [16], the new imposed constraints imply that finding experimentally viable 
areas, with a sizable cr^o_ pi becomes nearly impossible. Thus, for the present study 
we will focus on case (i). We will go through the same examples as in the previous 
subsection. 

Let us start with the regime where M\ m [/,, and as an example ensuring compatibil- 
ity with WMAP, choose A x = 400 GeV, A K = -200 GeV, fi = 130 GeV, with tan/3 = 5 
(corresponding to what was already depicted on the right panel of Fig. 3 and the left 
of Fig. 7). As shown in the previous subsection, there exist regions in the parameter 
space where the neutralino fulfils all experimental constraints and has the correct relic 
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density. The latter are characterised by neutralinos with a significant singlino fraction 
and/or a small mass. In this case, one of the allowed regions is close to the tachyonic 
area and exhibits very light singlet-like Higgses, potentially leading to large detection 
cross sections. This is indeed the evidenced on the left-hand side of Fig. 9, 

where the theoretical predictions for cr^o_ p are plotted versus the lightest neutralino 
mass. The resulting <x^o_ p spans several orders of magnitude, but, remarkably, areas 
with cr^o-p ~ 10 -7 pb are found. These correspond to the above mentioned regions 
of the parameter space with very light singlet-like Higgses (25 GeV < m h o < 50 GeV 
with Sf 3 > 0.99). The neutralino is a mixed singlino-Higgsino state with Nf 5 ks 0.35 
and a mass around 75 GeV. 

The sensitivities of present and projected dark matter experiments are also depicted 
for comparison. The small area bounded by dotted lines is allowed by the DAMA 
experiment in the simple case of an isothermal spherical halo model [17]. The larger 
area also bounded by dotted lines represents the DAMA region when uncertainties to 
this simple model are taken into account [18]. For the other experiments in the figure 
only the spherical halo has been considered in their analyses. In particular, the (upper) 
areas bounded by solid lines are excluded by EDELWEISS [21] 6 and CDMS Soudan 
[20]. Finally, the dashed lines represent the sensitivities of the projected CDMS Soudan 
and XENON IT [23] experiments. 

On the right-hand side of Fig. 9 we show the resulting cr^o_ p when the /x parameter 
and A K are varied to ji — 150 GeV, A K = 0, for which the effect of the different 
constraints on the (A, k) plane were represented on the right-hand side of Fig. 7. Since 
the areas of the parameter space with very light Higgses are ruled out by experimental 
constraints the detection cross section is not as large as in the previous examples. In 
the regions consistent with both experimental and astrophysical constraints the lightest 
Higgs mass is in the range 80 GeV < m h o < 120 GeV, thus leading to cr^)_ p < 5 x 
10~ 8 pb, within the sensitivity of projected dark matter experiments, such as CDMS 
Soudan. 

Let us now investigate the effect of changing the neutralino composition by mod- 
ifying the bino mass. As commented in [16], the largest values of the neutralino de- 
tection cross section were obtained for a mixed singlino-Higgsino composition, when 
/i < Mi < M 2 . In order to enhance the Higgsino composition we will consider examples 
where Mi is increased with respect to the /i-parameter. Such neutralinos annihilate 

6 Since the exclusion area due to ZEPLIN I [22] is similar to EDELWEISS we have not depicted it 
here, nor in any subsequent plot. 
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Figure 10: The same as in Fig. 9 but for an example with M 1 = 330 GeV, tan/3 = 5, 
A\ = 570 GeV, A K = -60 GeV, with fi = 160 GeV (left) and M 1 = 500 GeV, tan/5 = 5 
A x = 400 GeV, A K = -150 GeV, with fi = 130 GeV (right). In both examples, the 
resulting a^ USY is outside the experimental 2a region. 

more efficiently, thus leading to a reduced Q^o h 2 , so that the astrophysical constraint 
becomes more stringent. Nevertheless, as seen in the previous subsection, it is still 
possible to find areas of the parameter space with the correct relic density while simul- 
taneously fulfilling all experimental constraints. These regions corresponded to light 
singlet-like Higgses, which can potentially lead to sizable detection cross sections. 

First, Mi = 330 GeV will be taken, for an example with fi = 160 GeV, A x = 570 
GeV, A K = —60 GeV, and tan (3 = 5. The parameter space for this case was represented 
in Fig. 8, were we showed the effect of resonant annihilation channels on the allowed 
regions. The theoretical predictions for neutralino direct detection are shown in Fig. 10. 
In this plot, the various resonances appear as funnels in the predicted cr^o_ p ^ or ^he 
regions with the correct Q^o h 2 at the corresponding values of the neutralino mass 
(m^o ~ Mz/2 and m^o ~ m h o/2). Below the resonance with the Z boson, light 
neutralinos are obtained m^o < Mz/2 with a large singlino composition which have 
the correct relic abundance. The lightest Higgs is also singlet-like and very light, 
leading to a very large detection cross section, cr^o_ p > 10~ 6 pb. This corresponds to 
the allowed area of the (A, k) plane which lies in the vicinity of the tachyonic region in 
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Figure 11: Scatter plot of the neutralino-nucleon cross section as a function of the neutralino 
mass (left) and as a function of the lightest CP-even Higgs mass (right) for an example with 
tan/5 = 5, and the remaining parameters in the ranges 0.01 < A, k < 0.7, 110 GeV < 
M 2 < 430 GeV, -300 GeV < A K < 300 GeV, -800 GeV < A x < 800 GeV, and 
110 GeV < /i < 300 GeV. All the points represented are in agreement with LEP/Tevatron, 
ct^ USY , and BR(b — > S7) constraints, and have a relic density in agreement with the 
astrophysical bound (grey dots) or the WMAP constraint (black dots). 

Fig. 8. 

Remember however that these two examples with a larger bino mass were dis- 
favoured by the resulting muon anomalous magnetic moment, as it was illustrated in 
Fig. 1. 

One more example, this time for M x = 500 GeV, pi = 130 GeV, A x = 400 GeV, 
A K = —150 GeV, and tan f3 — 5 is represented in Fig. 10 and shows how large detection 
cross sections can also be achieved for heavier neutralinos. In this case (whose param- 
eter space was illustrated and discussed in Fig. 8) neutralino detection cross sections 
as large as cr^o_ p m 10~ 5 pb are possible while fulfilling experimental and astrophysical 
constraints. Once more, the occurrence of light singlet-like Higgses is crucial for en- 
hancing (J^o_ p and the sizable singlino component of the lightest neutralino (A^ 5 0.9) 
reduces the annihilation cross section and ensures the correct relic density. 

Finally, we show in Fig. 11a scatter plot of the theoretical predictions for the 
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neutralino-nucleon cross section as a function of the neutralino mass and the lightest 
Higgs mass when A, k, Mi, the \x parameter and the trilinear terms, A\ and A K are 
varied while keeping tan (3 — 5. In order to satisfy the a^ USY constraint, a small slepton 
mass, rriL,E = 150 GeV, has been used. Only the points which are in agreement with 
LEP/Tevatron, BR(b — > sj), and a^ USY limits and which, in addition, are consistent 
with the astrophysical bound or the WMAP constraint on the relic density are plotted. 
We clearly see how large detection cross sections are correlated to the presence of very 
light Higgses (m h o < 50 GeV). Neutralinos fulfilling all constraints and within the reach 
of dark matter experiments are possible with 50 GeV < m^o < 130 GeV. The upper 
bound on the neutralino mass is due to the lightest stau becoming the LSP. If the 
slepton mass is increased, heavier neutralinos can be found but the resulting af SY is 
soon outside the experimentally allowed range. 



5 Conclusions 

We have extended the systematic analysis started in [16] of the low-energy parameter 
space of the Next-to-Minimal Supersymmetric Standard Model (NMSSM), studying 
the implications of experimental and astrophysical constraints on the direct detection 
of neutralino dark matter. We have computed the theoretical predictions for the scalar 
neutralino-proton cross section, cr^o_ p , and compared it with the sensitivities of present 
and projected dark matter experiments. In the computation we have taken into ac- 
count all available experimental bounds from LEP and Tevatron, including constraints 
coming from B and K physics, as well as the supersymmetric contribution to the muon 
anomalous magnetic moment, a^ USY . Finally, the relic abundance of neutralinos has 
also been computed and consistency with astrophysical constraints imposed. 

We have found very stringent constraints on the parameter space coming from low- 
energy observables. On the one hand, a^ u is generally very small unless very light 
slepton (m,L t E ^, 200 GeV) and gaugino masses (Mi < 210 GeV) are considered, and 
slepton trilinear couplings modified in order to increase the LR mixing in the smuon 
mass matrix. On the other hand, the contribution to BR(6 — > s'y) is sizable, mostly 
due to the smallness of the charged Higgs mass, so that regions with tan (5 < 3 are 
disfavoured. 

Regarding the neutralino relic density, regions of the parameter space can be found 
where fl^o h 2 is in agreement with the WMAP constraint. This is possible when either 
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the neutralino mass is small enough for some annihilation channels to be kinematically 
forbidden or when the singlino component of the lightest neutralino is large enough to 
suppress neutralino annihilation. 

Remarkably, some of the regions fulfilling all the experimental and astrophysical 
constraints display very light Higgses, m h o ~ 50 GeV, which are singlet-like, Sf 3 > 0.9, 
thus allowing a sizable increase of the neutralino-nucleon cross section. Neutralinos 
with a detection cross-section within the reach of dark matter experiments are therefore 
possible, and have a mass in the range 50 GeV < m^o < 130 GeV. These neutralinos 
have a mixed singlino-Higgsino composition and are therefore characteristic of the 
NMSSM. 
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